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Abstract 

We introduce a "nonlinear" version of the Hamiltonian formalism 
which allows a self-consistent description of theories with degenerate La- 
grangian. A generalization of the Legendre transform to the case when 
the Hessian is zero is done using the mixed (envelope/general) solutions 
of the multidimensional Clairaut equation. The corresponding system of 
equations of motion is equivalent to the Lagrange equations, but contains 
"nondynamical" momenta and unresolved velocities. This system is re- 
duced to the physical phase space and presented in the Hamiltonian form 
by introducing a new (non-Lie) bracket. 
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Our previous result [T] was in generalizing the Legendre transformation to 
singular (degenerate) Lagrangians (with zero Hessian matrix) . For that a mixed 
(general/envelope) solution of the multidimensional Clairaut equation was in- 
troduced. 

In this paper we apply the above idea to construct a self-consistent "nonlin- 
ear" version of the canonical (Hamiltonian) formalism and present an algorithm 
to describe any singular Lagrangian system without introducing constraints [2] ■ 
To simplify matters we use coordinates, but all the statements can be readily 
converted to coordinate free setting [3l |4] . We also consider systems with fi- 
nite number of degrees of freedom. This is sufficient to explore the main ideas 
and constructions (this can be rendered to a field theory, e.g., using De Witt's 
notation [5]). 

First, recall the standard Legendre (-Fenchel) transform for the theory with 
nonsingular Lagrangian [6] . We then show its relation to the Clairaut equation 

in some details [1], which will be used to explain the main idea below. Lei0 
L A = 1, . . .n, be a Lagrangian given by a smooth functiorH of 2n 

variables (n generalized coordinates q"^ and n velocities = = dq^ /dt) on 
the configuration space TM. By definition, a Hamiltonian H as a dual 

function to the Lagrangian (in the second set of variables pa) constructed by 
means of the Legendre(-Fenchel) transform has the forrrH 

H{q^,PA)^snv[pBv'' -L{q^,v^)], (1) 

where the supremum is taken with fixed g"* and pA- In doing the Legendre 
transform, the coordinates q"^ are treated as fixed (passive) parameters of the 
duality transformation, and velocities, are independent functions of time. 
Then ([!]) leads to the supremum condition 

dL{q\v^) 

To obtain a dual function H we need to get rid of dependence on 

velocity in the r.h.s. of ([1]). This can be done in two ways: 

1) Direct way: resolve the condition ^ directly and obtain its solution as 
a set of functions — [q^,pA), then substitute them to ([1]) and obtain the 
standard Hamiltonian on the physical phase space T*M (see e.g. (6j) 

H^' {q^,PA) - PbV^ {q^,PA) ~ L (q^, (q^,p^)) . (3) 

This can be done only in the case of convex Lagrangian function (in the second 
set of variables u"*), which is equivalent to the Hessian being non-zero 



det 



d'^L{q^,v^) 



+ 0. (4) 



^We use indices in arguments, because by type of index we will distinguish them below. 

■^We consider time-independent case for simplicity and conciseness. 

•^We use summation convention for indices which are not in arguments of functions. 
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2) Indirect way: differentiate botli sides of ([U by momenta and use the 
suprcmum condition ([2]) to obtain the "dual supremum condition" in the form 

v-^V-{,-,p,)^'J^if^. (5) 

Then we substitute these velocities to (jU, which results in no manifest depen- 
dence of . Thus we obtain a partial differential equation with respect to 
Hamiltonian which in fact is the multidimensional Clairaut equation 

We call the transformation defined by ([5]) a Clairaut duality transform (or the 
Clairaut-Legendre transform) and i?*^' {<1^,Pa) a Clairaut-Hamilton function. 
Note that ^ is normally treated as a definition of dynamical momenta pA, 
but we should distinguish them from the parameters of the Clairaut duality 
transform pA'- before applying the supremum condition ([2]) they are assumed 
noncoincidental. 

The difference between the above two approaches is crucial for singular La- 
grangian theories. We thus label the resulting Hamiltonians by different indices. 
Specifically, the Clairaut equation ([5]) has solutions even in the case when the 
Hessian (j4]) is zero. So the Clairaut duality transform is more general and in- 
cludes the ordinary duality (Lcgendre-Fcnchel) transform as a particular case. 
To show this and find solutions of the Clairaut equation (O, we differentiate it 
by Pa to obtain 



PB - 



d^H-^ {q^,PA) 

' ~ dpA 



0. (7) 



So we have two possibilities depending on which multiplier in ([7]) is zero: 

1) Envelope solutions defined by the first multiplier in ([7]) being zero, this 
demand coincides with the supremum condition So we obtain the standard 
Hamiltonian ^ 

HtL [i^.Pa) \pa=va = H^' {q^,PA) ■ (8) 

2) A general solution defined the "dual Hessian" being zero 



d'H^^{q^,PA) 



dpAdpB 



dH-'{q^,PA) J, 



This gives — — c and then the general solution acquires the form 

OPA 

Hfen {q^^Pa) - PBC^ - L {q^,c^) , (10) 

where are arbitrary smooth functions considered in the Clairaut equation 
([6]) as parameters. Note that i?gg„ (g'^,^^) is always linear in the variables pA 
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which now do not coincide with the dynamical momenta, because we do not 
have the supremum condition ([2]). 

Now consider a singular Lagrangian theory for which the Hessian ([¥]) is zero. 

This means that the rank of Hessian matrix Wab = g^Ag^B \s r < n, and 
we suppose that r is constant. We rearrange indices of Wab in such a way 
that a nonsingular minor of rank r appears in the upper left corner. Represent 
the index A as follows: \i A — 1, . . . ,r, we replace A with i (the "regular" 
index), and, if yl = r + 1, . . . , rt we replace A with a (the "degenerate" index). 
Obviously, AeiWij ^ 0, and rankWij — r. Thus any set of variables labelled 
by a single index splits as a disjoint union of two subsets. We call those subsets 
regular (having Latin indices) and degenerate (having Greek indices). 

The standard Legendre transform is not applicable in the singular case be- 
cause the condition (|4|) is not valid [4j. Therefore the supremum condition ([2]) 
cannot be resolved under degenerate A, but it can be resolved under regular 
A only, because det Wij ^ 0. On the contrary, the Clairaut duality trans- 
form given by (O independent of the Hessian being zero or not Ij. Thus, 
we state the main assumption of the formalism we present here: the ordinary 
duality of convex functions can he generalized to the Clairaut duality for func- 
tions with zero Hessian. This can be rephrased by saying that the standard 
Legendre (-Fenchel) transform of nonsingular Lagrangian theory is generalized 
to the Clairaut-Legendre transform, and in both cases the corresponding trans- 
formation is described by the same Clairaut equation ©. 

To find its solutions, we again differentiate ^ by pA and present the sum 
([T]) in i? in two terms: regular and degenerate ones 



d^H-^ {q^,PA) 

dpAdpi 



dL {q^,v' 

Pa ~ 



dpAdpa 



As det Wij ^ 0, we suggest to replace (fTTj) by the conditions 



= 0. 
(11) 



f,- = 0- 13 

OpAOPa 

In this way we obtain a "mixed" envelope/general solution of the Clairaut 
equation, which can be also treated as a "partial" Legendre transform [1] . 
After resolving of under regular velocities = y ((7'^,pi, w") and writ- 

dH'^^lqA^pA 

ing down a solution of (|13p as — = (where are arbitrary 

OPa 

functions, unresolved velocities) we obtain a "mixed" Clairaut-Hamilton func- 
tion 

KL {q^,P^,Pa,v'') ^P^V' {q^,p^,v") +P^v"-L{q^,V' {q^,p.,v") ,v'^) , 

(14) 
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which is the desired Clairaut-Legendre transform written in coordinates. Note 
that ([14]) coincides with the "slow and careful Legendre transformation" of [8] 
and with the "generalized Legendre transformation" of [S]- 

1 he standard Lagrange equations oi motion ^ — ^"q;^ — - — — ^^g^ — - m our 
notation have the form 



'dt ^ 



dh„{q^,p,) _ dL{q^,v^) 



dq^ 



dt 



dq'^ 



where 



(q^,Pi 



dL (q 



(15) 



(16) 



The functions ha {q^,Pi) are independent of the unresolved velocities v" since 
rank VF^B = r. One should also take into account that now ^ = V'' {q"^,pi, w") 
. Note that before imposing the Lagrange equations (fT5| the 



and ^ 



arguments of L {q^^ v^^ were treated as independent variables. 

A passage to Hamiltonian formalism can be done by the standard procedure: 
consider the full differential of both sides of (|14p and use the supremum condition 
(fT^. which gives (till now the Lagrange equations of motion were not used) 



dH. 



I 

mix 



dp, 
dH, 



= V' {q' 



I 

mix 



dpa 



dm 



dq 
dH. 



I 

mix 



dq° 



,Pi,v 



dL iq' 



dL iq' 



dq" 



+ hp iq 



A M dv^ 



, f A s-1 dv^ 
V(9 'POJ q^- 



An application of the Lagrange equations (fT5|) yields the system of equations 
which gives a Clairaut-Hamiltonian description of a singular theory 



dH^ix 


dq^ 


dpi 


~ 'dt' 


dH"^ 

^ rmx 


dq"" 


dpa 


~ "dT' 


dH^ix 


dpi 


dq'- 


dt 


dH^ix 


dha {q^ 


dq"" 


dt 



+ hp [q- 



+ [pp + hp {q^,p^)] 



dv^ 
dq"' 



(17) 
(18) 
(19) 
(20) 



6 



This system has two disadvantages: 1) It contains the "nondynamical" mo- 
menta pa] 2) It has derivatives of unresolved velocities w". To get rid of them, 
we introduce a "physical" Hamiltonian 

Hq {q^,p,) = H-J^,^ {q^,p,,p^,v'') - [pp + h0{q^,p,)]v^ (21) 

Using (fT2|) and (|14p . one can show that the r.h.s. of ((2T|) indeed does not depend 
on "nondynamical" momenta pa and unresolved velocities . Then from (|17p - 
((20|) we obtain the system of (first order differential) equations which describes 
a singular Lagrangian theory 

^ = {q\Ho {q^,p^)} + {q\ h, {q^,p.)}v'', (22) 

^ = {p,,Ho {q^,p,)} + {p^.hp (g^,p.)} (23) 
Fap{q'^,p^)vf^ ^DaHo{q^,p^), (24) 

where {X, F} — f^rfl- — f|r|^ is the "regular" Poisson bracket (in regular 
variables). We introduce a "g"-total derivative" 

dX 

and a "g"-non-Abelian field strength (curvature)" 

r, ( A ^ a/tg (g^,P») dhp{q'^,pi) ^ \ U fA M /0R^ 

Fapiq ,P^)^ + {ha{q ,P^),h|3{q ,p.)|. (26) 

The system is equivalent to the Lagrange equations of motion due 

to our construction. 

In the case rankF^^ (^q^,pi^ = n — r, all the velocities v" can be found from 
(|24)) in a purely algebraic way. If rankF^^j {q^,Pi) — rp < n — r, then a singular 
theory has n — r — rp gauge degrees of freedom. In the first case one can resolve 
(Pi)) as follows 

^DaHo{q^,P,)F''^ {q^,p,), (27) 

where F"'' {q^,Pi) is the inverse matrix to Fq/3 {q^,Pi) , i.e. Fafi {q^,Pi) F^'^ {q^,Pi) = 
(52. Substitute ((27)) to ((22)) ~ ([23l) to present the system of equations for a singular 
Lagrangian theory in the Hamiltonian form as follows 

^^{q^,Ho{q^p.)}^, (28) 



dt 

^ = {p,,Ho{q^,p.)},, (29) 
where we define a new bracket 

{X,Y}p = {X,Y} + {X,ha (9^,^.0}^"" ('Z'^'PO D^Y. (30) 
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Note that the time evolution of any function X of dynamical variables 
{<l^:Pi) is also determined by the bracket ((30)) as follows 



The bracket pop is not anticommutative and does not satisfy Jacobi iden- 
tity. Therefore, the standard quantization scheme is not applicable here. We 
expect that some more intricate further assumptions should be made to quantize 
consistently singular systems within the suggested approach. 

To conclude, we describe Hamiltonian evolution of singular systems us- 
ing n — r + 1 functions of dynamical variables Hq Pi) and ha {q^TPi), a 
nonlinear Hamiltonian formalism. This is done by means of the generalized 
Clairaut-Legendre transform, that is by solving the corresponding multidimen- 
sional Clairaut equation. All variables are set as regular or degenerate. We 
consider the restricted phase space formed by the regular momenta pi only. 

There are two reasons why degenerate momenta pa are not worthwhile to 
be considered in a singular Lagrangian theory: 

1) the mathematical reason: there is no possibility to resolve the degenerate 
velocities as can be done for the regular velocities in (fT2)) : 

2) the physical reason: momentum is a "measure of movement" , but in 
degenerate directions there is no dynamics, hence — no reason to introduce the 
corresponding momenta at all. 

Thus there is no notion of constraint [21[T^ as restriction on "nondynamical" 
momenta, because eventually we do not consider the latter — thus nothing to 
constrain. Under this approach, the degenerate coordinates work as parame- 
ters analogous to n — r time variables (with n — r corresponding "Hamiltonians" 
ha {q^,pi), see psp). The Hamiltonian form of the equations of motion ([28]) - 
(|^ is achieved by introducing a new bracket ([SO]) depending on the above 
n — r+l functions. This bracket is responsible for the time evolution. However, 
is not anticommutative and does not satisfy Jacobi identity, and therefore its 
quantization requires non-Lie algebra methods. 

Thus, we presented here a general construction. Further details and exam- 
ples will appear elsewhere. 




(31) 
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